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We provide a realization of a singularity-free inflationary universe in the form of a simple cosmo- 
logical model dominated at early times by a single minimally coupled scalar fleld with a physically 
based potential. The universe starts asymptotically from an initial Einstein static state, which may 
be large enough to avoid the quantum gravity regime. It enters an expanding phase that leads 
to inflation followed by reheating and a standard hot Big Bang evolution. We discuss the basic 
characteristics of this Emergent model and show that none is at odds with current observations. 
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I. INTRODUCTION 

The idea that the universe we inhabit might be in 
a state of eternal inflation is not new. In fact it was 
realized quite soon after its initial proposal that inflation 
is usually future eternal; that is, in most cases there will 
always be regions of spacetime that are inflating in the 
future. So a natural question arises as to whether the 
universe was always inflating in the past; could inflation 
be past eternal also? Recent work by several authors 
01 seems to indicate that the answer to this intriguing 
question is decidedly "no"! Using Penrose-Hawking- 
Geroch techniques, it is argued in detail that a spacetime 
that (i) is past causally simple, {ii) is open, (iii) is 
described by Einstein's equations with a matter source 
that obeys the weak energy condition, and (iv) allows 
for inflation to be future eternal, cannot be past null 
complete. Crucial to the argument is the assumption 
that the universe is open or at least flat. As is often 
pointed out though, we are emerging into an era of 
"precision cosmology" and measurements of temperature 
anisotropics in the Microwave Background are able to 
place greater constraints on the curvature of the universe 
than ever before 0- The recent WMAP data seems to 
point to a universe that is close to {but not quite) flat, 
with a total density parameter of ^Itot = 1-02 ± 0.02. 
Nevertheless, until we have something like a three 
sigma signal, it remains premature to conclude with any 
certainty whether we live in a flat universe or not |34| . 
In particular, the WMAP data does not rule out closed 
models. In this light, it would seem that we have a 
promising avenue around the arguments of . If we live 
in a universe that is closed (albeit only marginally so) 
today, then it was always closed, and perhaps inflation 
is past-eternal after all. 

Recently arguments were put forward for several infla- 
tionary cosmologies that were past-eternal while avoiding 
any quantum gravity regime Each of the proposed 
models is spatially closed and described only with general 
relativity, ordinary matter, and minimally coupled scalar 
fields. Their existence argues against the suggestions of 
01 that inflationary universes are necessarily bounded 



in the past, and shows that a quantum gravity domi- 
nated era for the universe may not be inevitable; some- 
thing that has also been noted in The Eddington- 
Lemaitre cosmology is a well-known example of a uni- 
verse that is not past geodesic-incomplete, because of 
its positively curved spatial sections. Harrison has also 
given an exact solution with similar properties 0. His 
is a radiation-dominated closed universe with a positive 
cosmological constant. It starts from an Einstein static 
state, with a radius determined by the value of A, before 
entering a never-ending period of de Sitter expansion. 
However both these models do not exit inflation. Here, 
following 3] , we consider a universe filled with a dynam- 
ical scalar field, which is past asymptotic to an Einstein 
static model with a radius determined by the field's ki- 
netic energy. This model enters a period of de Sitter 
inflation that comes naturally to an end as the scalar 
field starts oscillating around the minimum of the poten- 
tial, before entering the standard hot Big-Bang expan- 
sion phase. Thus these are singularity-free inflationary 
universe modelsjlElj by-passing the restrictions of the sin- 
gularity theorems mentioned above |8fi| . They are finely 
tuned in terms of the initial conditions, although one can 
use entropy arguments to favor an initial Einstein Static 
phase for our universe |27l |. We consider possible impli- 
cations of this fine-tuning in the Conclusion. 



II. A POTENTIAL EMERGENT POTENTIAL 

Slightly more than a decade ago it was shown that 
by an inversion of the conventional viewpoint (begin- 
ning with a scalar field whose self interaction is dictated 
by some underlying particle physics considerations and 
subsequently determining the evolution of the universe) 
scalar field dynamics could be explicitly accounted for 
without a slow-roll approximation. Indeed a scalar field 
potential V{<j>) could quite easily be 'reverse engineered' 
for almost any desired behavior of the scale factor a{t). 
Several examples were explicitly computed and are sum- 
marized: 



2 



a{t) 




Aexp{ujt) 
A sinh(ajt) 
A cosh(wt) 


3^-1^2 + sinh^ f 2w((/« - M/B) 
3/4^1^2 + 52 sin2 (2uj{(t) - (l)o)/B^ 
(3n - 1)^2 exp ( ± 2(0 - 0o)/s) /2 



where k — SnG. The above respectively correspond to 
de Sitter exponential expansion, de Sitter expansion from 
a singularity, de Sitter expansion without a singularity, 
and power-law expansion respectively The eternal 
emergent universe, a nonsingular model past asymptotic 
to an Einstein static universe with topology R x S*"^ and 
radius oq 3> Lpi, can be treated in this way too. For a 
prescribed scale factor behavior of a{t) ~ oq + exp{ht) 
an associated potential is not too difficult to compute 
(see the appendix for details). While it certainly pro- 
duces the desired early time behavior, the potential does 
not have a definitely zero minimum. This is an unde- 
sirable feature as it points to a rather large cosmologi- 
cal constant, although there are ways around it. Also, 
because of the limits on the integration in the reverse- 
engineering construction, a graceful exit from the infla- 
tionary regime is not always guaranteed. Consequently 
we need to find a universe with a potential that is essen- 
tially the same as that in the eternal emergent universe 
at very early times but then goes to zero at some finite 
value of the field. That is, we need to look for a poten- 
tial V{(f>) that matches onto the reconstructed potential 
(Fig. 0) as t — > —00, with a long flat plateau, but then 
has a vanishing minimum value. 

Fortunately we need not look too far. Inflationary 
models based on higher derivative curvature terms go 
back to the remarkably prescient work of Starobinsky [gl 
in which the de Sitter phase was driven by the trace 
anomaly of the energy momentum tensor (see also 0). 
Among the variants on the original Starobinsky model, 
i?2-inflation based on a Lagrangian of the form £ ~ 
R + aR^ exhibits a particularly elegant implementation 
of a de Sitter phase with a linearly decaying Hubble pa- 
rameter, H. The i?2 term in this action is effectively 
an additional scalar degree of freedom which may be ab- 
sorbed by the introduction of a (non-dynamical) scalar 
field [13. [13 Einstein gravity is restored by an ap- 
propriate conformal transformation but at the expense 
of a d/unamical scalar field with an interesting poten- 
tial [HI IT^ [T3 . [T5I I. Since we build our model from 
a related potential, it is worth seeing how this works. 
Starting from an i?2_]-|2odified action 



and 



R = 



1 

f]2 



(2) 



Substituting this into ^ gives 



S = 



XV -9 



7^ an 



6 J d*xy^n\/f,(^g'"'dM) 



(3) 



The last term may be evaluated by noting that the di- 
vergence V^X^ = {l/yf^)^fJ,{y/^X^'). With this. 



d x^-g—g^^'^df.ndM 



after a boundary term is discarded. The action for the 
_R2 model in the Einstein frame is then written as 
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i4 / — ~( ~ 6q;2 

^"^l^-(TTW 



r^d^Rd.R - 



(4) 



To make contact with a canonical form for scalar field 
actions llS|j, it is usual to rewrite the scalar degree of 
freedom offered by the scalar curvature as 1^9 := -\/31n(l-|- 
2aR) so that 



S 



'l{R-\g'"'d^^d,V- ^{e-^'^^ - 



4a 



The effective potential in is just the reflection of that 
in Fig. n about = 0. Note that the parameters de- 
scribing the potential are fairly rigidly constrained (with 
respect to a) by the conformal transformation j^S] . 

Since it's conception the R^ model has received a sig- 
nificant amount of attention, largely as a result of the fact 
that the scalar field driving the de Sitter phase arises so 
naturally and is not inserted "by hand" solely to provide 
the inflationary dynamics. Indeed it was shown in [T^ 
that the model quite naturally supports a transient pe- 
riod of inflation followed by a FRW universe. Constraints 
on the coupling constant a are imposed by requiring that 
density perturbations be of an appropriate magnitude. 
Consequently, 10^2^/^ 2 < q, < 1916^-2^ This in turn 
determines the height of the plateau of the potential to 
be of the order IQ-^^ -IQ-^"^ x M^^i This form of the 
potential will form the cornerstone of our construction of 
the Emergent universe. 



S 



R + aR' 



= / d'^x^/^ {l + 2aR)R-aR^ 



(1) 



deflne := 1 + 2aR and make the conformal trans- 
formation g^jy I— > 'g^i, = fl^g^i, so that \f—g — i~l'^y/—g 



III. DETERMINATION OF THE PARAMETERS 
IN THE POTENTIAL 

Emergent type universes can be realized by relaxing 
the constraints on the R^ effective potential and consid- 
ering a spacetime filled with a minimally coupled, single 
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FIG. 1: The Emergent potential for the parameter choice 
B = 1 is identical up to an overall rescaling to the efltective 
potential for the scalar field in the Einstein frame of the R^- 
driven inflation model after relabeling (j) as —ip. 

scalar field cj) with a potential of tlie general form 

V^V{(l>) = {Ae^^-C)^+D, (5) 

where A, B, C and D are constants to be determined by 
the specific properties of the emergent universe. Then, 

!/'((/.)= 2AB (yl e^"^ - C) e^"^ , (6) 

and 

V"i(j)) ^2AB^ {2Ae^'t' -C)e^'l' , (7) 

where a prime indicates a derivative with respect to 0. 
Therefore, the above potential has a minimum at 0o — 
(l/B) ln{C/A) with Vq = = 0o) = D. Consequently, 
if we want to set the minimum of the potential at the 
origin of the axes, we must choose A — C and D = 0. 
Note that zero minimum for V{(l)) guarantees that there 
is no residual cosmological constant. Then, expression 
© reduces to 

V{cf>)^A{e^^~l)\ (8) 

By definition, the emergent universe corresponds to a 
past-asymptotic Einstein-static (ES) model Q. This 
means that 

1^(0^ -oo)-^, (9) 

where gq is the radius of the initial static model [28l |. 
It is then clear that for oq S> Lpi, where Lpi is the 
Planck length, the model can avoid the quantum regime. 
To determine an additional parameter recall that an ES 
universe filled with a single scalar field satisfies the condi- 
tion = 2/KaQ = (/), where k — SnG 1^. Hence, for 
an initially ES state we require that A — 2/Ka§, which 
brings the expression of the potential down to 

n</>) = A(e^'-l)'- (10) 



FIG. 2: A schematic illustration of the scalar field evolution in 
the Emergent universe. After leaving its initial static state (1) 
the model enters a slow-rolling regime (see case (ia)), or it goes 
through an intermediate pre-slow-roU phase (see case(ib)) - 
(2). In either case the scale factor grows sufficiently quickly 
to mitigate neglecting the curvature effects. A period of slow- 
roll (3) infiation is followed by a re-heating phase (4)and then 
by the standard hot Big Bang evolution. 

Therefore, the basic properties of the emergent universe 
have already fixed three of the four parameters in the 
original potential given by Eq. Q. The remaining pa- 
rameter (B) will be determined by considering other as- 
pects of the model and in particular by looking into 
the density perturbation spectrum (see Sec. VIII below). 
Nevertheless, we can still determine, the sign of B by 
demanding, without loss of generality, that 

y ^ i:^ (e^^ _ 1) e^^ < (11) 

for -oo < < 0. This means that e^"^ - 1 < and 
consequently that _B > 0. Incidentally, on comparison 
with the corresponding potential of the action, this 
is seen to correspond to a choice of negative coupling 
constant a as required to avoid manifesting tachyons and 
singular perturbative behavior in the model |l5j |. 

IV. LEAVING THE EINSTEIN-STATIC REGIME 

As the universe leaves the ES state the evolution of the 
scalar field is determined by the Klein-Gordon equation 

+ 3i/0 + I/' = , (12) 

where H ^ a/a \s the Hubble parameter. Therefore, 
given that > 0, we have expansion only if <f) < —V. 
Since the potential drops with increasing (j), we may 
consider the following two alternative cases: 

(i) (j) < 0, namely a decelerated scalar field, where the 
friction for the deceleration comes from the expansion. 
In this case one expects the kinetic energy of (f> drops. 
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Then, since 0^ ~ V initially and 1/' ~ 0, a slow-rolling 
period with seems likely. 

(ii) < 4> < — V, that is a "slowly" accelerating scalar 
field. Then, one expects the kinetic energy of the scalar 
field to increase. Under these conditions, a slow-rolling 
regime seems unlikely. Therefore, we are left with case 
(i), which splits further into two subcases: 

(ia) V < (j) < 0, namely a "weakly" decelerated scalar 
field, monitored by the familiar expression associated 
with slow-roll inflation 

3H,p + V' = 0, (13) 

since \(f>\ < \V'\; and 
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FIG. 3: The initial evolution of the scale factor in the Emer- 
gent universe. 



(ib) (j) < V' < 0, that is a "strongly" decelerated scalar 
field, described by the following form of the Klein-Gordon 
equation 



+ 3H(f> = , 



(14) 



since \4>\ > \V'\. 

So, as the universe leaves the ES state, it may evolve 
in a number of ways depending on the relation between 
(j) and V'. Cases (i) and (ii) lead to expanding models, of 
which (i) is the alternative corresponding to the Emer- 
gent Universe. Of the two possible subcases of (i), the 
first leads immediately to the standard slow-rolling infla- 
tionary regime. In the next sections we will concentrate 
primarily on this particular case. Before we proceed fur- 
ther, however, we should make the following comment 
with regard to case (ib). It involves a strongly deceler- 
ated, that is a very slow-rolling scalar field. Nevertheless, 
Eq. (|14ll also implies that (p oc a^^. The latter suggests 
that, as (f) drops rapidly, subcase (ib) could also lead to 
the familiar slow-rolling inflation. We may test this pos- 
sibility by considering the reverse engineered potential 
given in the Appendix (see Eqs. (|A13|) . (jA14|) '). where 
(j) dominates over V just like in (ib). Then, the point 
where the two potentials deviate should determine the 
initial conditions for the Einstein field equations that de- 
scribe subcase (ib). 



THE DURATION OF THE SLOW-ROLL 
REGIME 



We solve numerically and plot the results for the scale 
factor in Fig. O The graph clearly shows that a has the 
familiar exponential increase associated with standard 
slow-rolling inflationary models while at earlier times ap- 
proaches a constant non-zero value. 

The duration of the slow-rolling inflationary regime, in 
a spatially flat model, is determined by the usual slow- 



roll parameters 0] 



V" 



and ^ 

(15) 

where Mpi is the Planck mass. Throughout slow-roll 
approximation, the above satisfy the constraints e{4>) <C 
1 and |?7(0)| <C 1, which provide the limits of the slow-roll 
regime. Applied to the potential of the emergent universe 
(given by Eq. with B > Q and for -oo < (/> < 0), 

the above constraints read 



and 



1^(0)1 



(eB0 - 1)^ 



2^2 (2c^^ - 1) 



(cB0 _ ly 



< 1 



(16) 



(17) 



respectively, where we have set Mpi = 1 for simplicity. 
It should be emphasised that, although our model is spa- 
tially closed, the effect of the curvature (which is dom- 
inant at very early times) becomes negligible after few 
e-foldings, and only re-emerges in the recent universe. 

After a rather lengthy, but fairly straightforward anal- 
ysis, one can show that neither constraint provides a 
lower bound for cj). Therefore, the slow-roll regime starts 
at an arbitrarily small value of (j). In practice, this means 
that the emergent model starts slow rolling at a few e- 
foldings after leaving its initial ES state, at a finite 0;, 
when the curvature effects have become negligible. Given 
that the static regime corresponds to (/) — s- — oo, is very 
small and the available number of e-foldings can be very 
large. 

Note that, by employing ifTCIl and ((T7|l . we can also 
show that for any positive value of B > there is always 
a negative value for (f> at which the slow-rolling regime 
ends i.e., where e, ~ 1. For example, when _B = 1 we 
find that < - ln(l + \/2) (for Mpi = 1). 



VI. THE NUMBER OF E-FOLDINGS 

For standard slow-roll inflation the number of e- 
foldings is given by the expression 



HAt, 



(18) 



where 0j, (pf and t;, tf are the initial and final values 
of <f) and t respectively, and H is the Hubble parameter. 
The latter is generally given by 0, 



<p + V 



(19) 



where fc = 0, ±1 is the 3-curvature index. After few e- 
foldings, however, the effect of the curvature term be- 
comes negligible, while for slow-roll inflation we have 
(j>^ V. Thus, throughout the slow- rolling regime 
H — ^ kV jZ. Combining expression (|18|l with this re- 
sult we obtain 



TV 



0/ y 



(20) 



given that At = d(/)/0 and that = — y/voAd^, as 
the slow-rolling version of the Klein-Gordon equation 
{i.e. with ~ 0) guarantees. 

Applying the above to the case of the emergent uni- 
verse, namely inserting the potential (|10|l . we arrive at 
the expression 



2B 



(21) 



which provides the number of e-foldings associated with 
the emergent universe as a function of the yet undeter- 
mined parameter B. This is easily computed numerically, 
and in Fig.Qlwe plot A^ against the value of (/) at the start 
of the slowroU period for various values of the parameter 
B. Clearly, depending on the value of 0i, sufficient e- folds 
are easily obtainable in this model. Note that A^, while 
certainly very large, is nevertheless finite since the value 
of (pi can be stretched back only to the point where the 
curvature effects become appreciable. From a different 
point of view, the reason there are only a finite number 
of efoldings is because the initial value of the scale-factor 
is non-zero (unlike the standard infiationary model with 
K = 0). 



PQ 




FIG. 4: The number of e- folds obtained during a slowroll 
regime in the Emergent model plotted against (pi and for 
varoius values of the parameter B. 



At the same time, the Klein-Gordon equation gives (p — 
— V'/VSkV, having set ~ 0. Combining the two we 
arrive at 



-da 



(23) 



Applied to the emergent universe, by using expressions 
(|10|l and Hll|) . the above lead to the differential equation 



1 At l-e-' 

a 2B~^ 



which integrated gives 



ln(^ 

, at 



(24) 



0-^(<^-<^0- (25) 



So, the value of the scale factor at the end of the slow- 
rolling regime depends crucially on (pi, namely on the 
value of (p at the onset of slow-roll infiation. As expected, 
the smaller pi is the larger the final value of a. 



VIII. THE DENSITY SPECTRUM 

For structure formation purposes it is crucial to de- 
termine the density contrast at horizon crossin g 5 e- 
foldings before the end of inflation. Following [l^ we 
have 



VII. SCALE FACTOR EVOLUTION 



6p 



Ho 



N=50 



(26) 



N=50 



Starting from the Friedmann equation of the slow- 
rolling inflationary regime (see Eq. H19|) l we have 



During the slow-rolling regime the Hubble parameter is 
given by Eq. H19() and the above takes the form 
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Substituting the potential of the emergent universe (see 
Eq. (|10|) '). we obtain the following expression for the den- 
sity contrast associated with the model 



6p' 

P / Hor 



ao 



1 



(28) 



Ar=50 



which also depends on the yet undetermined parameter 
B. In view of the COBE observations the above is con- 
straint by 



ao 



(eB0 „ if 



10" 



(29) 



N=50 



Consequently, the CMB anisotropy limits allow us to 
express B as a function of oq, the radius of the initial 
ES state, which therefore becomes the key parameter of 
the emergent model. In Fig. we plot Sp/p against ao 
for several values of the parameter B. The appropriate 
value of the Einstein radius is then read off from the in- 
tersection of the perturbation curves with Sp/p = 10^^. 
For example, the B = 1 curve satisfies the density 
perturbation requirement when ao — 10^ Lpi. Again a 
comparison with the corresponding potential shows 
that a ^ Qq ^ lO^^Mp;^, within the required range for 
that parameter This result fixes the remaining 

parameter in our initial potential (see Eq. ((SJ) leaving 
us with just one more free parameter; the scalar field 
value when slow-roll commences 

Having used the amplitude of fluctuations to constrain 
the parameters of the model, it remains to check that 
other quantities whose values are tightly constrained by 
CMB data remain within their experimental bounds. 
Among these we count the power and polarization spec- 
tra and the spectral index of the density fluctuations. 
As a preliminary check, we note that the spectral index 
of scalar perturbations may quite easily be estimated by 
following the prescription of [23|. We thus define the 
horizon-flow parameters ei and €2 through the usual slow- 
roll parameters e and 77 as 



ei e/87r 

62 := (e-27;)/4^ 



(30) 
(31) 



in terms of which the spectral index ng — 1 = 
d\nVk/d\nk = —2€i — €2- Using the expressions lfTB|l 
and (|17|) for the slowroU parameters, an estimate for the 
first horizon-flow parameter is obtained as ei ~ 0.00014. 
Subsequent computation of the spectral index yields 
|ns - 1| 0.012 or ns 0.987. This result, although 
only a rough estimate is certainly compatible with those 
of |3, ■ We expect that a computation of the power 
spectrum will not deviate significantly from the standard 
result either and hope to report on this in forthcoming 
work 

With the above value for the Einstein static radius, we 
anticipate that the CMB anisotropy spectrum will have 



Q. 

to 
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FIG. 5: Plots of 5p/p against the size of the Einstein static 
universe ao for various values of the parameter B. As B 
increases from B = 0.1 through 5 = 1 the value of ao for 
which Sp/p ^ 10"'' increases from ao ~ 6.02 x lO'^Lp; to 
ao ~ 0.42 X lO'^Lpi 



the usual Sachs- Wolfe plateau and peaks as confirmed 
by current observations Moreover, once the scalar 

field potential is fixed we have a residual one-parameter 
family of models - parameterized by 0i - all asymptotic 
to the same Einstein static universe and within which is 
contained a class of models with fltot = 1 + £ where e is 
arbitrarily small. Since these models are all closed and 
begin from the same state they serve well to illustrate 
the emergent model. However, the fact that ^tot is not 
forced to be 1.02 avoids both excessive fine-tuning of the 
Einstein static radius and a possible cosmic coincidence 
problem |40| . 



IX. AVAILABLE ENERGY FOR REHEATING 

Inflation starts at ¥{(/) —00) = 2/KaQ and ends at 
(j) = with Vo = 0, which means that the maximum 
energy "stored" in the inflaton field is 



(32) 



In other words, the maximum energy available for reheat- 
ing is inversely proportional to the square of the radius 
of the initial ES state. If Oq = lO^Lp; we find 



¥{(/) -> c») ~ lO^^Gey 



(33) 



Therefore, provided that thermal equilibrium has been 
achieved and that the reheating process is efficient, the 
temperature at the beginning of the standard Big Bang 
evolution can be as high as 



T 



(34) 



Obviously, for less efficient reheating T^^ < 10 °GeV 
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X. DISCUSSION 

Closed inflationary models have not had a happy time 
of late. Indeed many authors have taken the latest mea- 
surements of Microwave Background isotropies, in par- 
ticular the WMAP best fit estimate of Qtot = 1.02 ± 0.02 
to signal a flat infinite universe. Yet there are two ways 
to interpret these observations. It could be that the uni- 
verse is in fact flat and more precise experiments in the 
future will eventually whittle away at the error bars un- 
til ^tot = 1 to within, say, one part in a million. On the 
other hand, a comparison of BOOMERANG data ^ 
(with a best fit of ritot = 1.02t°;°^) with WMAP data 
would seem to indicate a convergence on fltot = 1.02 > 1 
with increasing resolution. Admittedly, the data set from 
which we draw this conclusion is rather limited and it is 
with much eagerness that we await results from the Euro- 
pean Space Agency's PLANCK satellite. However, until 
this debate is settled one way or the other, one is forced 
to take closed models seriously. In stark contrast to some 
claims in recent literature |2j, it is not too difficult to 
construct consistent, single field inflationary models in 
a closed universe. We give one such construction here 
(but invite the reader to see [2^ for another, simple yet 
remarkably elegant model). 

The "Emergent Universe" proposed in ;3I] is a simple 
closed inflationary model in which the universe emerges 
from an Einstein static state with radius gq 3> Lpi, 
inflates and is then subsumed into a hot Big Bang 
era. The attractiveness of the proposed model is that 
one can avoid an initial quantum-gravity stage if the 
Einstein Static radius is larger than the Planck length. 
One might then ask whether such a model has a simple 
representation and whether it lies within the boundaries 
of current observations. In this paper, we provide a flrst 
explicit construction of such a universe. As such, it is 
a manifestly nonsingular closed inflationary cosmology 
that begins from a meta-stable Einstein static state and 
decays into a de Sitter phase and subsequently into 
standard hot Big Bang evolution. Inspired by an exact 
reverse-engineered potential constructed as in (see 
the appendix), this phenomenological model employs a 
single scalar field with a potential very similar to that 
arising in conformally transformed i?^-inflation only 
with a relaxing of some of the rigidity of that potential. 
In particular, beginning with a four-parameter potential 
Q, we can immediately fix two of the parameters by 
fixing the origin of the potential and requiring vanishing 
cosmological constant. By requiring this potential to 
match onto the exact reverse-engineered one at early 
times {\(t>\ 3> 1), a third parameter of the potential is 
related to the radius of the initial Einstein universe by 
— 2/KaQ. The remaining two-parameter model is 
then shown to exhibit all the desired properties of the 
Emergent universe model. Among others, these include 
a sufficient number of e-folds to solve the late-time 
fiatness problem, a spectrum of density fiuctuations with 
magnitude of the order of 10^^ and sufficient energy in 



the scalar field to allow for adequate post-inflationary 
reheating. Of these we find that the spectrum of density 
fiuctuations provides another constraint that allows us to 
relate the remaining two parameters so that effectively, 
the only parameter in the model is the size of the 
Einstein static universe. As for the re-heating energy, 
we show that for an initial radius of ao = lO^Lp; we can 
obtain re-heating temperatures of around 10^^ GeV . 

Of course, an immediate question arises as to a viable 
mechanism that realizes the initial Einstein static uni- 
verse that our model emerges from. At this point several 
avenues present themselves. Among the more promising, 
we note two; the first of these is the observation of 
"29] that the Einstein static universe is one of only two 
asymptotic solutions of the Ramond-Ramond sector 
superstring cosmology field equations The size of 

the positively curved Einstein universe in this picture 
is controlled by the level number of the Kac-Moody 
algebra of the conformal fields living in the compactified 
internal space. Indeed, this observation leads quite 
naturally to the tantalizing possibility of realizing 
an Emergent-like universe within a string cosmology 
context [331. A second, equally intriguing, possibility is 
that the initial Einstein static universe is created from 
"nothing" by some quantum tunneling process |3ll l32j | . 
Indeed, finiteness of the tunneling action requires that 
the universe created through instanonic tunneling be 
closed 33j. It is not implausible, then, that through 
spontaneous quantum fluctuations, a closed universe 
could be created in a long lived but transient Einstein 
static state which then makes a transition to a finite 
lifetime de-Sitter and subsequent marginally closed FRW 
phase along the lines described above. It remains to be 
seen that such a mechanism can be concretely realized 
in any quantum cosmological framework. 

The attractive aspects of the Einstein Static solution 
as a preferred initial state for our universe have been 
considered in the past. In fact, Gibbons has argued 
for the higher probability of an Einstein Static initial 
phase based on the model's maximal entropy [23|. Cru- 
cially, once the universe finds itself near the static state it 
could remain there for an undetermined amount of time. 
This is guaranteed by the neutral stability of the Einstein 
Static model against inhomogeneous (either pure fluid or 
pure scalar field) perturbations 26, 27, 2^]. Typically, 
expansion away from the static solution will lead to in- 
flation followed by the standard Hot Big Bang evolution. 
However, the reader will probably by now have noted the 
large degree of fine-tuning that went into setting up the 
initial state from which the universe emerges. Indeed, 
the emergent model is a very special trajectory in the 
space of possible infiationary evolutions. We have shown 
existence of such models, but not that they are probable. 

Some may regard this as a deadly blow to these mod- 
els, but we believe the case is wide open. Firstly, we note 
that although the idea that the universe should be prob- 
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able (and so not fine-tuned) is the dominant paradigm 
in cosmology at present, there is no scientific proof that 
this has to be the case. This is an unproven and indeed 
unprovable philosophical assumption, which may or may 
not be true 4i\- It is equally conceivable that - as was 
taken for granted in the past - whatever process causes 
the universe to come into being prefers a high-symmetry 
state. It is then relevant that the Einstein-static model is 
the highest symmetry non-empty Robertson- Walker uni- 
verse, and so would be preferred by such a process. 

Secondly, the models presented here show one can 
avoid the initial singularity if initial conditions were fine 
tuned in the way remarked on above. We believe it likely 
that this is a generic result: that given the usual physics 
of inflationary fields in the early universe (i.e. avoiding 
the introduction of 'shadow matter' which violates the 
weak energy condition), there is either a singularity at 
the start of the universe '43l or a fine-tuned initial state. 
This may be the real philosophical choice facing us: to 
decide which is worse, a space-time singularity, with all 
that that entails, or a fine tuning of initial conditions. It 
certainly seems very difficult to (phenomenologically, at 
least) construct a model that avoids both, and it is useful 
to recall Wheeler's characterization of space time singu- 
larities caused by gravitational collapse as the worst crisis 
facing theoretical physics. Nowadays we do not perhaps 
take such singularities seriously enough. 

Models of the kind presented here are useful in terms 
of making clear the alternatives facing us: we can in- 
deed avoid both a singularity and the quantum gravity 
regime, without introducing any exotic physics; but there 
is a price to pay in terms of fine-tuning. From some 
philosophical standpoints the high symmetry of the ini- 
tial state may even be an advantage. 

APPENDIX A: AN EXACT SOLUTION 

Consider a Friedmann-Robertson- Walker (FRW) uni- 
verse containing a minimally coupled scalar field (f) with 
Lagrangian density 

C^lidcj^r-Vi<p)^\^'-V{cb) (Al) 

where the last equality follows from assuming spatial ho- 
mogeneity of the (j) field i.e. (j) = (j){t). With this assump- 
tion, the stress-energy tensor takes the form of a perfect 
fiuid with energy density and pressure 

P4> = \^^-V{<P) (A2) 

respectively 0|- The classical equation of motion for (f> 

that follows from variation of the action S = J d'^x 

is 

<^ + 3iJ^-f^ = 0, (A3) 



where the Hubble parameter H := a{t)/a{t). The Ray- 
chaudhuri field equation for the FRW model with scalar 
field matter source and its first integral, the Friedmann 
equation are 

3ij + 3i?2 = SttG {V{(t>) - 0^) (A4) 

and 

3i72 + 3A = 8^G(i./)2-f y(0)) (A5) 

respectively. These equations, together with H12|l for a 
closed dynamical system from which the evolution of the 
universe model is determined. However, it is important 
to note that the Klein-Gordon equation is auxiliary in the 
sense that any solution of (jA4|) and ljA5|) with nonvanish- 
ing (f) will necessarily satisfy (|12ll so that the dynamical 
system in fact contains only two independent equations 
[43. These may be combined jT)] to give the (more con- 
venient) equivalent set of equations 

With these equations at hand, the potential V{4)) is con- 
tructed by 

• Specifying the constant k and a particular (mono- 
tonic) function a{t) and computing the associated 
Hubble parameter H and H. 

• Checking that the constraint k/a^ — iJ > is met. 
This assures the positivity of (j)^ as is necessary for 
a neutral scalar field. 

• Specifying an initial condition (p^ for (/)(t) and inte- 
grating the second of eas. (|A6p to get <f){t). This is 
then inverted (where possible) to give t(<j)). 

• Substituting into the first of eas. ljA6|l to obtain 
V{t) = V{t{(l))) and subsequently 1/(0). 

Thus, as advertised, this algorithm associates to a 
specified scale factor a{t) a potential V[(j)) to give a 
model that exactly solves the classical infiationary field 
equations. At this point some comments are in order: 
First, the reconstruction of the potential is manifestly 
non-unique as should be clear from the sign ambiguity in 
the choice of root in the equation (jA6|) . Secondly, it is 
not clear how sensitive this algorithm is to the choice of 
initial conditions for <f> i.e. on where on the potential the 
scalar field resides when it begins to roll. While general 
shape should be independent it might be expected 
that the details of the potential (location of extrema 
etc.) would vary with different initial conditions. This, 
however, does not affect the argument. 
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FIG. 6: A parametric plot of the re-constructed potential cor- 
responding to the exact solution a(t) = ao + exp(/it) normal- 
ized with respect to the size of the Einstein static universe. 
This re-constructed potential exhibits the same shape as the 
effective potential of the B? model in the Einstein frame but 
is plagued by a rather shallow minimum (only about 91% of 
its asymptotic value.) 



The emergent model is essentially a modified version 
of the Eddington-Lemaitre universe - k = +1, past- 
asymptotically Einstein static, singularity-free, without 
particle horizons and ever-inflating (see Q for more de- 
tails). Before applying the potential reconstruction tech- 
nique to this model, a few general observations about the 
scenario beg attention. The Einstein static state contain- 
ing matter with energy density pi and pressure pi = w^pi 
(— l/3<it;i<l)is characterized by 



i(i-w,)p, + y((^) = 



1 



47rG' 



(A7) 



where Qq is the radius of the spatial sections of the 
universe. If the sole content of this universe is the scalar 
field, as is the case in this note, pi = and the second of 
ea. (jA7|) requires that the scalar field have non- vanishing 
(but constant) kinetic energy. In this model with no mat- 
ter initially present, the mechanism envisaged in has (j) 
rolling at a constant speed along a flat potential (V = Vi) 
from (p = — oo at t = — oo to (f> — a,t t — where the 
potential first rises (and inflation is initiated) and then 
drops to a minimum a,i 4> = (j)f where the value of the 
potential is Vf = A/SttG ^ Vi. The field is carried over 
the hill by its non-zero kinetic energy and slow rolls to- 
ward the minimum where its damped oscillations reheats 
the universe. What follows is an attempt to realise this 
scenario. 

To this end, and to facilitate numerical computations it 
will prove useful to rescale quantities of interest as 



t/M, 



PI 



m 

ait)- 



■ Mpi(j){t), 
a{t)/Mpi, 



(A8) 



where M|,; := SttG. In these units the Raychaudhuri 
equation and its first integral, the Fricdmann equation 
become 



3H + SH^ = y(0)-02. 



(A9) 



which, together with the Klein-Gordon equation form the 
dynamical system governing the evolution of the closed 
scalar field dominated FRW universe. The corresponding 
equations ljA6|l for the potential and (j>{t) that are the 
starting point for the reconstruction of the potential take 
the form 

V{c^{t)) = H + 3H^ + \ 

4>'{t) - 2(\-h). (AlO) 



Following [3] consider the scale factor 
a{t) = A + Bexp{ht) , 



(All) 



where A, B, h are all positive constants. This universe is 
past asymptotic to an Einstein static phase, since a{t) — s- 
A as t —00. Thus, A is identified with the radius ao of 
the Einstein static universe. At late times, on the other 
hand, a{t) — > i3exp(ft,t) and the model approaches a de 
Sitter expansion phase. The second of eqs. IjAlOp then 
gives 




(A12) 



(ao + bf ' 

where b := B exp(M). The integral is easily evaluated to 



give 



nao 



1 — \/l — h'^aob\ 
1 + Vl - h'^aob) 



(A13) 



and with the b parameterization the potential is written 
3{hb)^ + h'^aob + 2 



Vib) 



{ao + by 



(A14) 



At this stage, an expression for V{4i) is, in principle, ob- 
tained by inverting eq. (jAlSjl and substituting b{(j)) into 
ea. (|A14|l . However the above form of ip{b) precludes such 
a simple treatment. Nevertheless, the general shape of 
V{(j)) is easily determined by a parametric plot. This is 
given in Fig El 
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